SetGeodeticParametersSystem Subroutine

SUBROUTINE SetGeodeticParametersSystem &
!
(system, prime_meridian)

USE Units, ONLY: &
! Imported parameters:
degToRad

IMPLICIT NONE

!Arguments with intent (in):
!!longitude, with respect to Greenwich, of the prime meridian
REAL (KIND = float), INTENT (IN) :: prime_meridian ![deg]
                                                   

! Arguments with intent (inout):
TYPE (CRS), INTENT (INOUT) :: system

!------------end of declaration------------------------------------------------

!set Geodetic parameters value
system % param (GEO_a)  = system % ellipsoid % a
system % param (GEO_b)  = system % ellipsoid % b
system % param (GEO_invf)  = system % ellipsoid % inv_f
system % param (GEO_prime)  = prime_meridian * degToRad !conversion to radians

!set Transverse Mercator parameters description
system % description (GEO_a) = 'semi_major_axis'
system % description (GEO_b) = 'semi_minor_axis'
system % description (GEO_invf) = 'inverse_flattening'
system % description (GEO_prime) = 'prime_meridian_longitude'

END SUBROUTINE SetGeodeticParametersSystem